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ABSTRACT 

Plasmas in an accretion flow are heated by MHD turbulence generated 
through the magneto-rotational instability. The viscous stress driving the ac- 
cretion is intimately connected to the microscopic processes of turbulence dissi- 
pation. We show that, in a few well-observed black hole accretion systems, there 
is compelling observational evidence of efficient electron heating by turbulence or 
collective plasma effects in low accretion states, when Coulomb collisions are not 
efficient enough to establish a thermal equilibrium between electrons and ions at 
small radii. However, charged particles reach a thermal equilibrium with their 
kind much faster than with others through Coulomb collisions, a two-temperature 
accretion flow is expected. We consider a Keplerian accretion flow with a constant 
mass accretion rate in the pseudo-Newtonian gravitational potential and take 
into account the bremsstrahlung, synchrotron, and inverse Comptonization cool- 
ing processes. The critical mass accretion rate, below which the two-temperature 
solution may exist, is determined by the cooling processes and the coUisional en- 
ergy exchanges between electrons and ions and has very weak dependence on the 
collision-less heating of electrons by turbulence, which becomes more important 
at lower accretion rates. Collision-less heating of electrons by MHD turbulence 
can no longer be ignored in quantitative investigations of these systems. 

Subject headings: acceleration of particles — accretion, accretion disks — black 
hole physics — plasmas — radiation mechanisms: thermal — turbulence 



1. Introduction 

Black hole accretion is one of the most powerful energy sources in the universe. When 
the luminosity of the system is close to the Eddington luminosity, the accretion can be de- 
scribed by the classical Shakura-Sunyaev disk (1973), which produces a multi-color blackbody 
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radiation with the flux and temperature determined by the black hole mass and accretion 
rate. Observations of galactic black hole X-ray binaries strongly support such a scenario 
(Gierlinski & Done 2004). The less energetic non-thermal high energy emission component 
frequently observed has been attributed to hot magnetized coronas above the disk (Gierlinski 
et al. 1999; Zhang et al. 2000). 

The magneto-rotational instability (MRI) has been generally accepted as the basic 
mechanism producing the less well-understood turbulence viscosity that drives the accre- 
tion (Balbus & Hawley 1991). The dissipated gravitational energy is first converted into 
MHD turbulence, which then heats the accretion fiow through viscous and Ohmic dissipa- 
tions. In these optically thick slim disks, the Coulomb collision time scales are much shorter 
than other relevant time scales, electrons and ions reach a thermal equilibrium so that the 
observed emission can be well-described with blackbody spectra modified by the radiation 
transfer through the disk structure. The acceleration and/or heating of electrons by turbu- 
lence can at most be studied with observations of the non-thermal component emitted from 
the collision- less coronas (Dermer et al. 1996; Li & Miller 1997). 

Many black hole candidates are often observed in a super- or sub- Eddington emission 
state, and there are many distinct observational phenomena, such as relativistic outfiows and 
quasi-periodic oscillations, that have not been well understood though many models have 
been proposed (Fender et al. 2004; Rcmillard et al. 1999; Falcke et al. 2004; Camenzind 
2005). Both theoretical investigations and observations suggest that, below a critical mass 
accretion rate M^r, a two-temperature accretion fiow likely develop near the black hole due to 
the much higher radiation efficiency of electrons than ions and inefficient Coulomb coupling 
between them so that the local thermal equilibrium can not be established between the two 
(Shapiro et al. 1976; Rees et al. 1982; Zdziarski et al. 2002; Titarchuk & Fiorito 2004). 
Such a two-temperature flow is also expected due to the fact that charged particles reach a 
thermal equilibrium with their kind much faster than with others through Coulomb collisions 
(Spitzer 1962). This paper studies the electron heating processes in these low states. 

Besides energy exchanges through Coulomb coUisions, charged particles can also be en- 
ergized by plasma waves through coUision-less processes. These energization processes will 
directly affect the characteristics of the observed emission and therefore play crucial roles 
in our study of these systems in the low states. As pointed out by Bisnovatyi-Kogan and 
Lovelace (1997), the Ohmic heating of electrons through dissipation of MHD turbulence can 
be very efficient. Sharma et al. (2007) recently found that electrons might also be energized 
in the dynamo processes of magnetic field amplification through the MRI. However, the 
complexity of processes in a turbulent plasma has made the collision-less electron heating by 
collective plasma effects a difficult problem (Rees et al. 1982). In these theoretical investi- 
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gations, certain assumptions have to be made on the couphng between charged particles and 
the turbulent electromagnetic fields to derive some quantitative results, which are usually 
sensitive to the prior assumptions. 

The collision-less electron heating by turbulent magnetic fields is often ignored in most 
phenomenological models (Shapiro et al. 1976). In the advection dominated accretion fiow 
models, it is simply assumed that a small fraction of the energy dissipated through viscosity 
is converted into electrons (Esin et al. 1998; Yuan et al. 2006), suggesting that this is a 
trivial process. Alternatively, one-temperature models have been proposed for Sagittarius 
A*, the compact radio source associated with the low-luminosity supermassive black hole in 
the Galactic Center (Schodel et al. 2002; Ghez et al. 2005), assuming that electrons and 
protons are coupled by turbulence effectively (Melia et al. 2000, 2001). 

On the other hand, the collision-less electron heating processes may be constrained by 
observations of the relevant systems. Recent studies of fiares from Sagittarius A* indicate 
that electrons can be heated efficiently by MHD turbulence and the distribution of relativistic 
electrons under the infiuence of a turbulent magnetic field can be approximated as relativistic 
Maxwellian (Baganoff et al. 2001; Genzel et al. 2003; Gillessen et al. 2006; Bittner et al. 
2007). The same processes may also play important roles in the heating of electrons in 
the two-temperature accretion fiows. In this paper, we consider a Keplerian accretion disk 
model and the dominant cooling processes in a fully ionized magnetized plasma. We show 
that Coulomb collisions with ions can not heat electrons efficiently and an extra electron 
heating process is required to explain observations of Sagittarius A* and the galactic X-ray 
binary Cygnus X-1 (Wilms et al. 2006). This is in contradiction with what is suggested in 
some of the previous studies. Magnetic turbulence could play such a role, and we believe 
that the collision-less heating of electrons by MHD turbulence should not be ignored in any 
quantitative studies of these low accretion states. 

The critical mass accretion rate is determined by the cooling processes and collisional 
energy exchanges between electrons and ions. Although the collision-less turbulence heating 
can reduce M^r by increasing the radiation efficiency and cooling rate, we find that this 
effect is insignificant. More detailed investigations of the radiation spectrum and magneto- 
hydrodynamic processes are needed to study these plasma processes quantitatively. 

The basic equations for the accretion disk are given in § [2l In § [3l we discuss the electron 
heating processes and show that observations of Sagittarius A* in the millimeter and sub- 
millimeter range are difficult to explain without introducing efficient electron heating by 
turbulence. The coohng processes are studied in §1] and the model is applied to galactic 
X-ray binaries in § [51 In § |6l we draw conclusions and discuss the model limitations and 
possible improvements in the future. 
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2. Basic Equations for Keplerian Flows in the Pseudo-Newtonian Potential 

We consider a fully ionized hydrogen plasma. Then the gas pressure and the thermal 
energy density are given, respectively, by 

P = nkB(Tp + Te), (1) 
£ = nkBil.5Tp + aTe) , (2) 

where the gas density, the proton and electron temperatures are denoted by n, Tp and T^, 
respectively, ks is the Boltzmann constant, and a = x[3Ks{x) + Ki{x) — 4K2{x)] / 4K2{x) 
with X = nieC^ /ksTg, where mg and c denote the electron mass and the speed of light, 
respectively, and Ki refers to the ith order modified Bessel function. Ki(x) — 2*"^(i — l)!/a;* 
as X ^0 and Ki{x) (7r/2x)^/^ exp(-a:;)[l + {Af - 1)/Sx] as x ^ oo. 

In the pseudo-Newtonian gravitational potential (Paczyiisky & Wiita 1980), the poten- 
tial and Keplerian angular velocity are given, respectively, by 



<P = (3) 

r — rs 

1/2 

(4) 



K 



GM 



r{r — r^Y 



where the gravitational constant, the black hole mass and the radius are denoted by G, M 
and r respectively, and rs = IGMjc? is the Schwarzschild radius of the black hole. 

Due to the decoupling between electrons and ions, the radiation efficiencies of two- 
temperature accretion fiows will not be as high as the optically thick slim disks. The accretion 
processes likely drive strong winds from the hot disk (Blandford and Begelman 1999), and 
there are currently no strong observational constraints on these processes. To simplify our 
model we consider the steady state properties of the accretion flow averaged in the vertical 
direction and assume a radius independent accretion rate, i.e., ignoring the effects of winds: 

M — —A'KrHvrn{'mp + rUe) , (5) 

where v,,. and H are the radial velocity and the scale height of the accretion flow, respectively, 
and nip is the proton mass. The vertical structure of the disk can be very complicated as 
suggested by MHD simulations (Hawley & Balbus 2002) and may be essential to explain the 
observed non-thermal high energy emission in low states (Dermer et al. 1996; Li & Miller 
1997). It, however, is not expected to introduce significant changes to our quantitative 
results below, which deal with the energetically dominant thermal emission component from 
the main body of the hot accretion fiow (Shapiro et al. 1976). 
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Following Blandford and Begelman (1999), we assume that the angular momentum flux 
through the disk is negligible, which is also expected for strongly magnetized disks (Agol & 
Krolik 2000). MHD simulations also show that the azimuthal velocity of the flow is given 
by the Keplerian velocity fx = rQi<i (Hawley & Balbus 2002), which requires the inner 
boundary radius > 3 rs, the radius of the last stable orbit. From the angular momentum 
conservation of the accretion flow, we then have 

^ dr \2r r — rs J 

where the kinematic viscosity 

_ 2/5p/?,A;B(Tp + r,)(r-rs)V/2 



(mp + me)(3r-rs)(G'M)V2 



(7) 



and in accord with Melia et al. (2001) f3p = (i?^)/87rP and f3i, is defined as the ratio of the 
total stress to the magnetic field energy density {B^/8tt) Then the radial velocity can be 
rewritten as 

^ PuPpkBjTp + Te){r - rs) , . 

~ (mp + me)(GMr)i/2 " ^ > 

Since the magnetic field is dominated by the toroidal component, one has 

>fcB(Tp + Te)(l + 2/5p^T'/' 



H 



GM{mp + m. 
GMA 

4vr/3,/5,[fcB(Tp + T,)]3/2(l + 2(3py/^r{r - rg 



(r - rs) , (9) 



GMM{mp + m,y/^ 



The energy conservation equation is given by 

^ {HrVr[Pil + 2l3p)+£ + n{mp + m,)[(l)-0.5vl + 0.5v^.]]} = -A, (11) 



Hrdr 



where A is the radiative cooling rate, and we have ignored any energy fluxes carried away 
from the disk by winds and waves and taken into account the effects of magnetic fields 
properly (Melia et al. 2001). Note that the work done by the torque force is given by 
—d[HrVrn{mp + me)vl^]/dr. Then we have 

dr VrU 



Note that the viscous parameter in the classical disk model a = Pi,j3p/{1 + 2f3p). 
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where 

£ = kB[Te{a + 1 + 2(3p) + Tp(2.5 + 2(3p)] + {rup + me)[0 - 0.54 + O.S^;^] (13) 

can be considered as the energy of the accretion flow per proton and Me /nip gives the 
outwardly directed energy flux through the accretion disk (Blandford & Begelman 1999). 



3. Electron Heating 

Because the viscous stress driving the accretion is induced by the MRI (Balbus & Hawley 
1991), the gravitational energy dissipation flrst generates MHD turbulence, which then heats 
the gas as the turbulence cascades from large scales to small scales. A fraction of the 
turbulence energy may also be carried away from the disk by MHD waves propagating toward 
large radii. However, this process is not expected to dominate the dynamics of the accretion 
flow at least for the relative more powerful states, when the hard X-ray luminosity of X-ray 
binaries can be a significant fraction of the Eddington luminosity. Very high accretion rates 
will be needed to produce the observed X-ray power if waves carry most of the dissipated 
gravitational energy toward large radii. The turbulence cascade also needs to be suppressed 
dramatically to make the wave escape process dominant. Therefore, in the steady state, the 
viscous heating rate 

,d[-(j) + O.dvl] 

r = Vrn[m„ + me) : 

dr 

should be slightly greater than the turbulence cascade rate 



SttH ' 

where C2 ~ 1 is dimensionless and determined by the properties of the turbulence and 

"A;B(T, + Te)(l + 2/3p)T^/' 



cs 



rrip + rrif. 



is the speed of fast mode waves (Liu et al. 2006b). There is therefore an intimate connection 
between the viscous stress and turbulence dissipation. It is, however, not clear how this 
energy is distributed between electrons and protons in a collision-less plasma and what 
determines this energy partition, especially for hot plasmas in the two-temperature accretion 
flow of black holes. In the following, we will show how observations of low states of black 
hole accretion systems may be used to constrain the coUision-less electron heating rate by 
turbulence. 
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The disk structure can be obtained once one specifies the electron heating and coohng 
rates and the electron and proton temperatures at the outer boundary. Following Blandford 
and Eichler (1987), we have the electron heating time by turbulence (Liu et al. 2006b) 

r.^^, (14) 

where Ci is a dimensionless constant, and 

2c{x + l) 



x'^K2{x) exp(a;) 



is the mean electron speed. Note that {ve)/c — > {SkBTe/ nmeC^y^'^ as x — > oo and the electron 
heating becomes relatively more efficient at smaller radii if electrons become relativistic and 
(ve) approaches to c. Then we have 

1^^^ + - — Z — ' 15 

where 

37rmemp(i'e)^ 
"""""^ " 256neMnA ' 

is the electron-proton energy exchange time through Coulomb collisions, e is the elemental 
charge unit, and InA ~ 15 for most astrophysical situations (Spitzer 1962). 

There are several ways to estimate Ci. For a non-radiative accretion flow, A = 0. At 
large radii, we can ignore the kinetic energy associated with the radial motion and a — 1.5. 
If we ignore the Coulomb coupling term, we have the following solution: 

cs = (A:B(Te + Tp)(l + 2/3p)/(mp + me))'/' = [3(l + 2/5p)/(5 + 4/5p)]^/V, (16) 
H = [3(l + 2/3,)/(5 + 4/?,)]V2^, (17) 

Vr - -[3/3,/3p/(5 + 4/3p)]7;K, (18) 

n = [3/{5 + 4pp)]-^/\l + 2Pp)-^/^M/47iPaU^ip + meyvK, (19) 

k^Tg = [3/{5 + 4f3,){l + 2(3,)]-'[n/72C',f3X]my^, (20) 

VTe = [3/(5 + 4/?p)]2[72C2/32/52/vr(l + 2/5p)](mp + me)/me-l 

= 1.5 X lQ^[Ci(3M{l + 0.8/3j,)]V(l + 2/3p) - 1 . (21) 

Then we have 

Ci ^ 8.1xlO-^[{Tp/Te+l){l+2/3p)]'/^{l+0.8Pp)/P,Pp > 8.1x10-3(1+2/3^)^/2^1+0.8/?^) . 
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It is interesting to note that the proton temperature satisfies an equation similar to equation 

dr Tacfr ' 

suggesting that the heating time of protons is the same as electrons. Since the mean speed 
of protons (vp) is usually different from (ve), it is likely that the dimensionless constant Ci, 
which is determined by the microscopic physics of resonant and/or non-resonant coupling 
between particles and turbulence, is different for protons and electrons. Otherwise, Tp/T^ = 
nip /me, and we have 

Ci = [(5 + APp)/18P,PpMl + 2/3p)/2]i/2 . (22) 



For this solution, the Coulomb energy exchanges between electrons and protons (with 
Tcoui independent of the radius r) become relatively more important at larger radii. So we 
have a one-temperature flow at large radii, which develops into the above solution at the 
radius, where the Coulomb energy exchange time becomes comparable to the accretion time: 

^ _ (5 + 4Ppy/^ 



The transition radius is therefore given by 

7c^memp{mp + me)GMc^[l + 2(3pf\ ( Ti^mpll + 2(3pf L^ddV''^ 



rt/rs 



3[5 + 4/3p]2i5/2e4 in XCfp.ppM J + 0.8/3^] 2^3/2 me In XCfp.pp Mc^ J 

) " ^ '° I [1 + ) [m^j 



\[l + 0.8/3p]Cf/3,/3pMcy 

where i^Edd = 4:iiGM{mp + me)c/aT is the Eddington Luminosity and ax = (87r/3)e'^/mgC^ 
is the Thomson scattering cross-section, and we have used the above solution to obtain the 
electron temperature and assumed Te = Tp to get the last expression. We see that two- 
temperature flows can readily develop at small radii even without considering the electron 
cooling processes. 

When Coulomb collisions become less efficient at smaller radii, in principle the distri- 
butions of electrons and protons are not necessarily Maxwellian. However, at a given mean 
energy, the time scales for protons and electrons reaching equilibrium with themselves are. 



(24) 



^With the assumptions adopted above, this equation actually can be derived from equations (|lip and 
(|15p . It is therefore not an independent equation. 
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respectively, ~ {rne/mpY^'^TconX and ~ {rne/mp)Tconh which are much shorter than rcoui, the 
time scale for electrons and protons reaching thermal equilibrium with each other (Spitzer 
1962). The electron and proton distributions therefore may start to deviate from Maxwellian 
at {me/mpY/^rt and {me/mpY^^rt, respectively. So the transition radius, where the proton 
distribution may deviate from Maxwellian, is more than ten times smaller than the radius, 
where the two-temperature flow starts to develop. The transition radius for electrons is 
even smaller, and a Maxwellian distribution can be achieved at small radii through relativis- 
tic effects and cooling processes (Wolfe & Meha 2006; Bittner et al. 2007). We therefore 
expect that the electron distribution be always dominated by a thermal component and, 
over a significant radius range, a two-temperature fiow develop self-consistently. The fact 
that the low-hard state spectra of X-ray binaries can be fitted with the thermal inverse 
Comptonization model also suggests a dominant thermal electron component. 

At very low accretion rates, the Coulomb collision time scales can be much longer than 
other relevant time scales near the black hole, and we do not expect Maxwellian distribu- 
tions for both electrons and ions. This is a completely new regime, where the behavior of 
the turbulent plasma has not been well studied. In these cases, the "temperature" of the 
particles should be interpreted as a measurement of their mean energies. The actual electron 
distribution may be constrained by studying its radiation spectrum. We note that the proton 
distribution is always expected to deviate from Maxwellian at very small radii, where the 
energy dissipation rate is high, and a high energy proton population may be responsible for 
the formation of jets in the low states (Dermer et al. 1996; Liu et al. 2006a). A detailed 
investigation of these effects is beyond the scope of this paper. 

Prom r > and the above solution, we have Pi, > 2C2/S. The viscous stress therefore 
needs to be comparable to the turbulent magnetic field energy density. Since the viscous 
stress should not be much greater than the turbulent magnetic field energy density (Pessah et 
al. 2006), the fact that C2 ~ 1 requires that Pg can not be much less than P. So the energy 
carried away by winds and/or waves should account for a small fraction of the dissipated 
gravitational energy except that the turbulence cascade from large scales to small scales is 
suppressed for some reason so that C2 becomes much less than 1. 

We also note that for this solution, t^c oc r^/^. The turbulence electron heating effectively 
converts a constant fraction S of the viscously dissipated gravitational energy into the internal 
energy of electrons. The electron heating rate by turbulence is given by ankBTg/Tac- Then 

S ~ [7r(5 + 4/3p)(l + 2Pp)/216C^PlP^]mJ{mp + m,)<l 

gives CiPi^Pp > 6.3 x 10-^[(1 + 0.8/3p)(l + 2/?p)]V2. Although this result is similar to that 
proposed in the advection dominated accretion fiow models (Esin et al. 1998), the electron 
heating studied here becomes very efficient when cooling reduces the electron temperature or 
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the electrons become relativistic since the heating time is proportional to the mean velocity 
of electrons. Accordingly, Coulomb collisional heating may suppress the turbulence heating 
by increasing (fe). 

Since the heating is more efficient when electrons become relativistic, an even stricter 
constraint on Ci can be obtained by considering the accretion processes at small radii. For 
A = and Tg ^ Tp, a = 3, the accretion time is given by 

^ _ (8 + Ap,yii^ 

'''''''^ ?,Mp{GMyn- ^^^^ 

To avoid electrons being overheated, the electron heating time needs to be longer than the 
accretion time, then we have: 

^,,„,,>(i±i^)"'^(^)'%0.3, (26) 

where we have used the fact that, for the accretion flow we are considering, r > 3rg. 

For A = 0, we solve equations (1121) and (1151) numerically to constrain the heating rate 
more quantitatively. For (3p = 1, we find that CiP^Pp needs to be greater than 0.15 to 
avoid the proton temperature decreasing below zero at small radii at certain accretion rates 
due to efficient electron heating. Panel "a" of Figure [T] shows the temperature proffies for 
M = 10"'', 10~^, lOLEdd/c^- (Note that the Eddington accretion rate is usually defined 
as MEdd = lO^Edd/c^O The other model parameters Ci = 1.5, Pu = 0.1, Pp = 1, and 
Te = Tp = GMrrip/lOkBro at the outer boundary = lO'^r^. When M is small. Coulomb 
collision effects are negligible. We recover the solution given by equations (fT6] — [21!) at large 
radii. At small radii, the electron temperature goes beyond that given by equation (1201) 
due to the relativistic effect discussed above. With the increase of M, the electron and 
proton temperatures become closer to each other at larger radii. At small radii, because of 
efficient heating of relativistic electrons, the electron temperature is higher than the proton 
temperature for M = 10~^LEdd/c^- When the accretion rate becomes comparable to or 
higher than the Eddington accretion rate ~ lOLEdd/c^, the electron and proton temperatures 
become identical due to very efficient Coulomb collisions. 

The linearly polarized millimeter and sub- millimeter emission from Sagittarius A* (Aitken 
et al. 2000; Schodel et al. 2002; Chez et al. 2005) reveals a mass accretion rate below 
10~^LEdd/c^ (Quataert & Gruzinov 2000). X-ray observations and detailed modeling sug- 
gest an even lower accretion rate (Baganoff et al. 2001; Melia et al. 2000; Liu & Melia 2002). 
VLBI observations have shown that this emission comes from a region within ~ 15 rs of the 
black hole and should originate from a hot accretion torus (Shen 2005; Melia et al. 2001). 
The dotted line in Panel "a" of Figure [1] shows the electron (lower) and proton (upper) 
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temperature profiles for M = 10~^LEdd/c^ and Ci = 15, wliicli is ten times bigger tlian tfiat 
for tlie otlier lines, giving mucli less efficient electron heating by turbulence. Although we 
have adopted a relatively high accretion rate. Coulomb collisions are not efficient enough to 
heat electrons to relativistic energies required to produce the observed synchrotron emission. 
A fit to the millimeter and sub-millimeter spectrum and polarization gives Cij3yj3p ~ 0.1, 
which is comparable to the maximum heating rate discussed above, where cooling is absent, 
and is more than 10 times smaller than that given by equation (122|) (Liu et al. 2007). One 
can show that Ci for electrons is more than 10 times smaller than that for protons in this 
case. This could be due to the high mean momentum of protons so that the action of protons 
on turbulence in their stochastic scattering processes is important. The proton heating is 
then suppressed giving rise to a more than 10 times longer acceleration time than electrons 
at a given mean speed of both particles. 



4. Electron Cooling 

For a fully ionized magnetized plasma, synchrotron, inverse Comptonization (IC), and 
bremsstrahlung are the dominated emission mechanisms (Rybicki & Lightman 1979). In the 
low accretion states of black holes, the disk is optically thin to bremsstrahlung radiation. 
We therefore have the bremsstrahlung cooling rate 

Abrem = TTT ^T? 9^ = 1-4 X IQ-^'T^/'ti' , (27) 

where qb ^ 1.2 is the Gaunt factor and h is the Planck constant. For equations (fT6|) - (!2T|) . 
Abrem oc r"''/^ and the bremsstrahlung cooling time of electrons Tbc = a/cB^^e/Abrem oc r. 
Including this energy loss process, we find that the disk collapses to a low temperature flow 
toward small radii when M > = 2.63 x lO^^LEdd/c^ for an outer boundary radius of 
lO^rs. This can be understood by comparing Tbc with the accretion time r^is oc r^/^ [See eq. 
Then we have 



.^^^ ^ 3V3..,/.,../.gy;.K ( 3 y ^ 2.9(^3/.)"=^ , (28) 

where we have assumed that Tg = Tp, i.e. Ci = [(5+4/5p)/18/3i,/3p][7r(l+2/9p)me/(me+mp)]-^/^, 
and used the model parameters for Panel "a" of Figure [H to obtain the last expression. The 
critical mass accretion rate decreases with the increase of the outer boundary radius. 

Panel "b" of Figure [1] shows the temperature profiles for several values of M. The other 
model parameters are the same as Panel "a" . Note that the bremsstrahlung cooling has little 
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effect on the temperature profile for M <^ Mcr because, compared with the accretion time 
given by equation fl23|) . this coohng becomes less important toward smaller radii. For M = 
Mcr, we see that the temperatures are lowered by the cooling near the outer boundary. Since 
the cooling rate is proportional to the density square, the cooling time scale is proportional 
to M, and McrC^/Ledd is independent of M. Pair production will increase the cooling rate 
and decrease Mcr- These results are similar to those given by Narayan & Yi (1995) though 
electron heating by turbulence has been included in our calculations, implying that the 
collision-less electron heating does not affect the bremsstrahlung cooling rate at large radii 
significantly. 

One may define a critical radius with equation fl28l) : 



A hot accretion flow can only exist below r^r- When the transition radius rt given by equation 
f l24|) is greater than r^, as is the case for relatively high accretion rates, a hot two-temperature 
flow may develop below r^. For low accretion rates, can be much greater than rt, we have 
a hot one-temperature flow between and r^r, which develops into a hot two-temperature 
flow below Vf. The thin lines in Panel "b" corresponds to this case, where ~ lO^r^ and 
r„ = IOV5. From n = r^r, we have Mc^/L^aa ^ 156/32/32/(1 + 2/3^)^4(1 + 0.8/3^)^/2 

~ 117rs(l + 2/3p)^/2/(i + 0.8/3p). However, how exactly the hot flow may develop 
toward small radii is still an open question, which may depend on the nature of the large-scale 
flow and the coupling between the disk and its corona (Meyer- Hofmeister et al. 2005). 

When the source is optically thin to synchrotron radiation, the cooling rate is given by 

^syn = 7r^(P )B ^ 1.06 X 10 i^nfi2 , (30) 

where (^2) is the mean momentum square of the electrons (Melia & Coke 1999). For equations 
f lTB]) - fl^ and X <^ 1, Agyn oc r^^. The corresponding cooling time scale oc r^/2^ which is 
more important at smaller radii. For x 3> 1, Agyn oc and t^c oc r^/2_ However, most of 
the thermal synchrotron emission is emitted at 

x + 20 60(1 + 0.05x)e572 ^ „^ 

ue ~ -^—i^c = -^r^ TT^ = 8.4 X 10^^7,2 1 + 0.05x / X + 1 Hz , 

x + 1 47r(x + l)meC ''^ ' 

where 7c = l/x-|-l (Liu et al. 2006). The thermal synchrotron emission (emissivity per solid 
angle) and absorption coefficients are given, respectively, by (Petrosian 1981; Mahadevan, 
Narayan, & Yi 1996) 

x/3e^(p2)(a;+l) ^ 

= Q 3 4/ 2\n I n in -. BuzMlizM) , (31) 

87rm;fc4(7^)(l + 0.19x) 
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_ SuC^[exp (hu/kBTe) - 1] ^ ■Ken{p'^) I{zm){x + 1) , , 

~ 2/IZ/3 " 3v^7>e(72)SfcBTe;2M(l + 0.19x) ■ ^ 

where (7^) ~ (x + 12)(x + l)/x^ is the mean Lorentz factor square of electrons, (x + 1)/(1 + 
0.19x) is chosen to take into account the non-relativistic effects (see below), and 

I{zm) = 4:.0505zl^^\l + 0A0z~}^'^ + 0.53l6zM^^)exp{-1.8899zli^), (33) 
zm = z//z/c = 47rmeCz//3ei?7f . (34) 

We then have the optical depth through the emission region 

/ N _ f ^rr T^e{p^)n H I{zm){x+1) 

Tyiu) = Ky[l']H = — ; -, (35) 

^' ^' 3v^m,(72)fcBre7,^S^Af(l + 0.19x) ' ^ ' 

which, for a given zm-, is proportional to r^^/^ for equations f|T6|) - fl2T]) and x ^ 1 and 
proportional to r^^^ for x ^ 1. When Ty{iJE) becomes greater than unity at large radii, one 
has to take into account the self-absorption effects. The corresponding synchrotron cooling 
rate may be approximated as 

^ _ 16\/37r2A;BTeZ/| (72)(1 + x)(l + O.lQx) 

~ 452/[(x + 20)/(x + l)]c2iJ 72(1 + 0.05x)2 

SOx/Se^fe Te53 (72)74(1 + 0.05x)(l + O.lQx) 

37r/[(x + 20)/(x + l)]m3c5 i7(x + 1)2 

~ . . 1 n-i3^ei?' + 12)(x + 1)3(1 + 0.05x)(l + 0.19x)/(20) 
_ 5.4 xiu — a;6/[(x + 20)/(x+l)] ' ^"^^^ 

where the numerical factor is chosen so that the cooling rates given by equations fl30p and 
fl36|) are equal when t^^ue) = 1, and Agyn — SnkBTei'^/Sc'^H for x ^ 1. That is, for 
cyclotron radiation in the optically thick regime, we assume a black body spectrum cutting 
off at ve = t'c- This cooling is important at small radii since Agyn oc r~^^/^ and Tsc oc r^'^/^ 
for equations ffTB]) - fl^ and x <^ 1. For x ^ 1, A^yn oc r~23/4 g^^d r^c oc r^^/^. 

We will assume that the IC is in the Thomson limit and the seed photons for IC are 
provided by the synchrotron radiation. Then we have the total cooling rate 

A = Asyn + Aic + Abrcm = Asyn[l+87r(A-Abrcm)^/c52]+Abrem (37) 



A,y,(l-87rAsyn///cS2)-i+A 



brem 



Panel "c" of Figure [T] shows the temperature profiles with the above cooling processes in- 
cluded for M = 3.4 X 10^ Mq. The model parameters are the same as Panel "b" except 
M = 2.63 X 10-6 (thick), 2.63 x 10"^ (medium), 2.63 x 10-2 (thin) LEdd/c^- We note that 
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the electron temperature never goes above ~ 10^^ K, which may provide an alternative ex- 
planation to the observed upper limit in brightness temperature of powerful extragalactic 
radio sources (Readhead 1994). The critical mass accretion rate Mcr = 2.63 x 10~^LEdd/c^ is 
the same as that obtained above, where only the bremsstrahlung cooling is considered. This 
is consistent with the fact that the synchrotron and IC coolings are unimportant at large 
radii. At small radii, these cooling processes do reduce the electron temperature. 

Because the synchrotron cooling rate in the optically thick regime is proportional to 
n^/'^/H, the cooling time is proportional to M^/^. This cooling becomes relatively more 
efficient for smaller black holes so that the electron temperature decreases with the decrease of 
M for given Mc^ / L^dd- also decreases slightly with the decrease of M. For M = SAMq, 
we find Mcr = 2.13 x 10~^LEdd/c^- Panel "d" of Figure[I]is the same as Panel "c" except that 
M = 3.4M0. The fiatten part of the electron temperature in the inner region for M ~ M^t 
is dominated by Coulomb coUisional heating and optically thick synchrotron cooling, which 
is very sensitive to the electron temperature when x < 1. We therefore expect that, for 
M = Mcr, the intrinsic brightness temperatures for extragalactic radio sources are higher 
for bigger black holes and there should be a correlation between the brightness temperature 
and the source luminosity. Pair production effects need to be incorporated to give more 
quantitative predictions. 

Due to the low radiation efficiency at large radii, the critical mass accretion rate is not 
very sensitive to the electron heating rate either. Panel "e" of Figured] show the temperature 
profiles for Ci = 15. The other parameters are the same as Panel "d", and Mcr = 2.51 x 

10-2LEdd/c2. 

5. Application to X-ray Binaries in the Low-Hard States 

It has been suggested that Coulomb collisions are efficient enough to heat electrons to 
explain observations of galactic X-ray binaries in the low-hard states (Shapiro et al. 1976; 
Esin et al. 1998). Although some of the relevant models may achieve an electron temperature 
in line with observations by adjusting parameters describing the source structure and/or 
dissipation processes, the related assumptions are not well justified both theoretically and 
observationally. It is also not obvious whether these models can explain the recently observed 
anti-correlation between the electron temperature and the Thomson optical depth of Cygnus 
X-1 over a large dynamical range (Wilms et al. 2006). The model proposed here is self- 
consistent and has the basic parameters to describe the dynamics and radiative processes. 
It is therefore well positioned to uncover new physical processes when applied to specific 
observations. 
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With the parameters adopted above, the model clearly can not explain the high hard X- 
ray luminosities observed in some of the galactic X-ray binaries. The hard X-ray luminosity 
can be a significant fraction of L^dd for some sources. A much higher M^r is needed to produce 
such a high luminosity. To increase M^r, one may increase the viscosity by increasing f3iy or 
decrease the outer boundary radius. The virial temperature at r = IC^rs is already a few 
tens of keV, and we are interested in showing how electrons reach this temperature range 
from a large-scale cool accretion flow. We will consider the first option and set the outer 
boundary temperature /cs^e = k^Tp = 5 x lO^^nipC^ ~ 5 keV in the following. Panel "a" 
of Figure [2] shows the temperature profiles for M = 34:Mq, Ci = 0.2, Pp = 1, and = 3. 
The critical mass accretion rate M^r = 21.9LEdd/c^; which is consistent with equation fl28|) . 
Therefore, if the hot accretion fiow develops naturally with the processes considered here, a 
very high viscosity is needed to produce a hard X-ray luminosity on the order of O.lLEdd- 
This statement is not necessarily true if the formation of hot accretion fiows is triggered by 
some instability near black holes. 

However, for M ~ Mcr, the electron temperature is more than 1 MeV near the black 
hole, which is much higher than the observed values in the low-hard states (See Fig. [2] a). 
It is clear that, besides the local internal cooling processes considered above, other cooling 
processes have to be introduced to bring the electron temperature to the desired energy range. 
Observations of X-ray binaries in the low-hard states show that photons from the large-scale 
optically thick slim disk can provide sufficient soft photons to cool electrons in the inner 
two-temperature fiow through IC (Barrio et al. 2003). The cooling rate will also increase if 
we consider a global radiation transfer and include the cooling due to pair production. Here 
we treat these cooling processes approximately by assuming the energy density of the soft 
photons produced outside the two-temperature accretion fiow is three times higher than the 
local magnetic field energy density so that there is an extra cooling term 3 times higher than 
that given by equation fl30l) . Panel "b" of Figure [2] shows the temperature profiles with this 
cooling effect included. In this case, Mcr = 0.933LEdd/c^, which is clearly affected by this 
extra cooling process. The other model parameters are the same as Panel "a" . 

According to this model, the IC of the external soft photons dominates the cooling, 
and the electron temperature is nearly constant at small radii due to the balance between 
Coulomb collisional heating and this cooling for M = M^. The value of the electron tem- 
perature depends on the ratio of the external photon energy density to that of the magnetic 
field 77: 

, ^ / 9 In Xrrie \ 9 , . 

ksTe ~ —= — rrieC^ . (38) 

\4:\/2'Kj3primpJ 

The electron heating by turbulence becomes relatively more important for lower values of M. 
However, with M, which is proportional to the Thomson optical depth of the accretion fiow. 
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varying by two orders of magnitude, the electron temperature in the inner region remains 
about 30 keV. This is not consistent with the recently observed anti-correlation between the 
electron temperature and the Thomson optical depth of Cygnus X-1 in the low-hard and 
intermediate states (Wilms et al. 2006; Titarchuk 1994). 

By increasing the electron heating rate by one order of magnitude, i.e. with Ci = 0.02, 

we obtain temperature profiles in Panel "c" . Because IC cooling is already important at the 
outer boundary, M^r = 0.891LEdd/c^ decreases slightly due to more efficient electron heating, 
which results in more cooling of the accretion flow. The other model parameters arc the same 
as Panel "b" . Wc sec that the electron temperatures arc in the observed energy range. Panel 
"d'' shows profiles of the Thomson optical depth r = 2aTHn and e for the same models in 
Panel "c". We quahtatively recover the observed anti-correlation between ksTg and r. The 
radiation efficiency of the accretion flow is ~ 0.01 giving rise to a luminosity of ~ 10^5-38 
erg s~^, which is also in line with observations. We note that CiP,^Pp — 0.06 implying very 
efficient electron heating by turbulence. 



6. Conclusions and Discussion 

To explain the millimeter and sub-millimeter spectrum and polarization of Sagittarius 
A* and the recently observed anti-correlation between the electron temperature and the 
Thomson optical depth of Cygnus X-1 in the low-hard and intermediate states, we show 
that electrons need to be heated efficiently by turbulence plasma waves. Coulomb collisions 
are not effective enough to heat electrons to the temperature required by observations of 
Sagittarius A* that has a very low accretion rate. 

The critical mass accretion rate, below which a two-temperature solution may exist, is 
determined by the radiative cooling and Coulomb coUisional heating processes and is almost 
independent of the electron heating by turbulence. To reproduce the observed high luminos- 
ity of X-ray binaries in the low accretion states, a high viscosity is required to increase the 
critical mass accretion rate, and the radial velocity of the accretion flow can be comparable to 
the Keplerian velocity. We show that the electron cooling for galactic X-ray binaries needs to 
be dominated by IC of an external soft photon fleld originated presumably from a large-scale 
cold slim disk. Synchrotron, synchrotron self-Comptonization, and bremsstrahlung processes 
inside the hot accretion flows are not efficient enough to cool the electrons to the observed 
temperature range. Assuming that the energy density of these external photons is 3 times 
higher than the local magnetic field energy density, we find that the electron temperature 
at the critical accretion rate is determined by this energy density ratio. A lower external 
photon energy density gives rise to a higher electron temperature. This high external photo 
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energy density suggests that the radiation field may have important dynamical effects and 
may be partially responsible for the required high viscosity. 

Recent observations of Cygnus X-1 show that the electron temperature in the IC model 
(Titarchuk 1994) decreases by more than a factor of 10 when the Thomson optical depth, 
which is proportional to the accretion rate, increases by two orders of magnitude in the low- 
hard and intermediate states. We show that very efficient electron heating by turbulence can 
reproduce such an anti-correlation. Without this extra heating process (besides the Coulomb 
collisional heating), the electron temperature is almost independent of the accretion rate. It 
appears that such an anti-corrcIation may also be explained by adjusting the properties of 
the external photo field (Esin ct al. 1998). It is, however, not clear whether the observed 
dramatic changes in both the electron temperature and the Thomson optical depth can be 
reproduced in such a scenario. Detailed modeling of the external photo field is required to 
address this issue and to fit the broadband spectrum of Cygnus X-1 in the low-hard and 
intermediate states. 

We have shown that a two-temperature fiow can develop self-consistently at small radii 
for low accretion rates mostly due to the fact that electrons and protons reach a thermal 
equilibrium with themselves much faster than with each other through Coulomb collisions. 

However, at very low accretion rates, there is no guarantee that the proton distribution 
is Maxwellian although electrons may reach a quasi-thermal distribution due to balancing 
between the cooling and heating processes. Observations of the low-hard states of X-ray 
binaries also suggest that the electron distribution is dominated by a thermal component. 
When Coulomb collision time scales are much longer than other relevant time scales, the 
particle "temperature" should be interpreted as a measurement of the internal energy but 
not a thermal distribution. 

With the pseudo-Newtonian potential and assumption of Keplerian azimuthal velocities, 
processes below ~ 3rs can not be described accurately by our model. These processes can 
have important observational consequences (Titarchuk & Fiorito 2004) and may drive strong 
outfiows responsible for the observed radio emission (Heinz & Sunyaev 2003; Fender et al. 
2004; Falcke et al. 2004). By dropping these assumptions and considering the momentum 
conservation in the radial direction, one may build a more complete model to study these 
processes (Popham & Gammie 1998). 

In the model we have ignored the vertical structure of the disk, dynamical effects of 
the radiation field, winds and/or outfiows, outwardly directed waves, and pair productions, 
which can be important for X-ray binaries. These may introduce a factor of a few changes 
in our model parameters, especially in regards to the temperature profile of protons and 
the transition from an one-temperature fiow at large radii to the two-temperature inner 
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region, which may have important imphcations on the origin of quasi-periodic oscillations. 
But our conclusion of the presence of efficient collision-less electron heating by turbulence is 
not affected. Electrons can also be accelerated by MHD turbulence to very high energies in 
coronas of the disk and produce the observed non-thermal high energy spectral component 
(Gierlihski et al. 1999). These electrons carry much less energy than the thermal electron 
population considered in this paper and can be treated with appropriate corona models 
(Dermer et al. 1996; Li & Miller 1997). 

Our treatment of the electron heating by turbulence is different from that proposed 
in the advection dominated accretion flow models (Yuan et al. 2006), where the energy 
equation is separated into two equations, one for electrons and the other for ions, and it is 
assumed that a fixed fraction of the viscously dissipated energy is deposited into electrons. 
We point out that equation (ITTl) is the more appropriate energy equation one should adopt 
and there is no reasonable justification to split it into two equations except that the coupling 
between electrons and ions through plasma waves or turbulence can be ignored. Our results 
and the MRI clearly show the opposite, i.e. turbulence plays crucial roles in an accretion 
flow. It is also difficult to understand why a fixed fraction of the viscously dissipated energy 
goes into electrons. Electron heating by plasma waves should depend on the properties of 
the plasma and turbulence (Bisnovatyi-Kogan & Lovelace 1997; Quataert 1998). Equation 
f|T5|) is perhaps more appropriate to describe the temperature evolution of electrons. 

It is very challenging to study turbulence, one of the commonest natural phenomena. It 
originates from the non-linear nature of complex many body systems. In astrophysics, MHD 
turbulence is ubiquitous, and one of the most important effects of it is the energization of 
charged particles in collision-less plasmas. The microscopic details of this stochastic particle 
acceleration process are still a matter of debate. Even the common practice of treating 
MHD turbulence as a spectrum of plasma waves is not well justified, especially for strong 
turbulence. The results in this paper provide a means to study the energization of electrons 
and ions by turbulence, specifically the energy partition between the two particles, with 
observations. Our results suggest that the electron heating time by turbulence is at least 
10 times shorter than that of protons at a given mean particle speed. Detailed modeling 
of the radiation spectrum and MHD simulations of the accretion fiows will lead to more 
quantitative results. 

We note that the dimensionless parameter Ci, which characterizes the electron heating 
by turbulence in our model, decreases by a factor of > 50 from the supermassive black hole 
of Sagittarius A* to the stellar mass black hole in Cygnus X-1. This may be related to the 
dramatic increase of the turbulence viscosity described by (3^, which changes from 0.1 for 
Sagittarius A* to 3 for galactic X-ray binaries. Indeed, the product of the two changes by 
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less than a factor of 2. Because (5i, > 2C2/3 is proportional to the energy dissipation rate 
of MHD turbulence, the energy partition between electrons and ions is described by C1C2. 
Our results suggest that in the two-temperature accretion states of black holes, this partition 
weakly depends on the black hole mass and mass accretion rate. 

This work was funded in part under the auspices of the U.S. Dept. of Energy, and 
supported by its contract W-7405-ENG-36 to Los Alamos National Laboratory. Some ideas 
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Feng Yuan for their hospitalities and Dr. Xue-bing Wu and Lei Qian for useful discussions. 
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Fig. 1. — a: Temperature profiles of electrons (solid) and protons (dashed) for M — 10"^ 
(thick), 10~^ (medium), 10 (thin) L-^dd/c^- Cooling is not inchidcd to demonstrate the 
effects of electron heating by turbulence and Coulomb collisions. The electron and proton 
temperatures at the outer boundary are chosen so that A-'eTg = A-'eTp is equal to 1/10 of 
the gravitational binding energy of protons and e < 0. Ci = 1.5, /ij, = 0.1, and /Sp = 1. 
The dotted lines have Ci = 15 and M = lO^^LEdd/c^- Note that the proton temperature 
decreases sharply toward small radii for M — O.OOlLEdd/c^ and the electron and proton 
temperatures are identical for M — lOLEdd/c^. See text for details, b: Similar to "a" with 
the bremsstrahlung coohng included and M — 2.63 x 10~^ (thick), 2.63 x 10"^ (medium), 
2.63 X 10~^ (thin) I/Edd/c^- The disk collapses to a cold disk at small radii for M > M^r — 
2.63 X 10~^I/Edd/c^(?'/10'^rs)~^/^. The dotted lines have Ci = 15 and the same accretion rate 
as the thick lines, c: Same as "b" with all the cooling processes included and M = 3.4 x 
IO^Mq. We see that the synchrotron and IC coolings only affect the temperature profile of 
electrons at small radii, d: Same as "c" but for M = 3.4M0. Because the synchrotron cooling 
time in the optically thick region scales as M^/^, the electron temperature decreases with M 
and Mc, = 2.13 x IQ-^LEdd/c^ e: Same as "d" but for Ci = 15. M^r = 2.51 x IQ-^LEdd/c^ 



- 25 - 




0.5 1 1.5 2 2.5 3 3.5 4 




0.5 1 1-5 2 2.5 3 3.5 4 



Log,o(r/rs) 




0.5 1 1.5 2 2.5 3 3.5 




0.5 1 1.5 2 2.5 3 3.5 4 



Log„(r/rg) 



Fig. 2.— a: Same as "e" of Figure □ but for M = 34M0, Ci = 0.2, /5p = 1, = 3 
and M = 2.19 x lO'^ (thick), 0.219 (medium), 21.9 (thick) LEdd/c^ In this case, Mcr = 
21.9LEdd/c^- b: Same as "a" but with an external soft photo field, whose energy density 
is 3 times higher than the local magnetic field energy density. M = 9.33 x 10~^ (thick), 
9.33 X 10-2 (medium), 0.933 (thin) LEdd/c^ and M^r = 0.933LEdd/c^ c: Same as "b" but 
for Ci = 0.02, and M = 8.91 x lO'^ (thick), 8.91 x lO'^ (medium), 0.891 (thin) LEdd/c^. 
Mcr = 0.891LEdd/c2- d: Profiles of e (dashed) and r (solid) for the three models in "c". The 
thickness of the lines are the same for a given accretion rate. The luminosity of the disk for 
the thin lines is more than 9% of L^dd- 



